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CN ■ ABSTRACT 

CO I We present photometry of the nearby galaxy NGC 5128 (Centaurus A) observed with the 

. PACS and SPIRE instruments on board the Herschel Space Observatory, at 70, 160, 250, 350 

^SJ ■ and 500 fim, as well as new CO J = 3 — 2 observations taken with the HARP-B instrument on 

I the JCMT. Using a single component modified blackbody, we model the dust spectral energy 

■ distribution within the disk of the galaxy using all five Herschel wavebands, and find dust 

temperatures of ^ 30 K towards the centre of the disk and a smoothly decreasing trend to 

20 K with increasing radius. We find a total dust mass of (1.59 ± 0.05) x 10^ M©, and a 
total gas mass of (2.7 ± 0.2) x 10^ Mq. The average gas-to-dust mass ratio is 103 ± 8 but 
' we find an interesting increase in this ratio to approximately 275 toward the centre of Cen A. 

H . We discuss several possible physical processes that may be causing this effect, including dust 

■ sputtering, jet entrainment and systematic variables such as the Xco factor. Dust sputtering 

by X-rays originating in the AGN or the removal of dust by the jets are our most favoured 
explanations. 

Key words: galaxies: individual: Centaurus A - galaxies: ISM - infrared: ISM - submillime- 
tre: ISM. 



1 INTRODUCTION 

Elliptical galaxies were long thought to be gas and dust poor, 
but in the past few decades, several observational surveys have 
shown that a significant fraction of these galaxies contain a rich 



* Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA, 
f E-mail: parkintj@mcmasterca 



interstellar medium (ISM). Detections of cold dust (temperatures 
of ~ 30 K) in non-peculiar early type galaxies using data from 
IRAS co uld only be made i n about 12 percent of the sam- 
ple (e.g. ISregman et alj [l998h : with improvements in our abil- 
ity to detect the cold dust component of the ISM, w e are find- 
ing m ore evidence for dust within eai i y-types (e.g. 'Kmi et al.l 
2004 ; IXilouris et al.||2004l : iLeeuw et al.ll20oi : lYoung et al.ll2009l : 
iSmithetal.ll20I2t) . Other studies suggest there is a relationship 
between a detection of cold dust in these ga laxies and the like - 
lihood of a radio source at the core (e.g. IWalsh et al.1 [198^ . 
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Cold ga s has also been found in elliptical galaxies, both molecu- 



lar (i.e. Thronson et al. 1989: Wiklind et al. 1995; Welch & Saed 


2003; lYoune et al, 2011), and atomic fe.a. iKnaoD et al. 


1985; 


Breaman et al. 


19921; iHuchtmeiei 1994; Huchtmeier et al 


1995; 


Moraanti et al 


Il999l,l2006l; Welch et alj2010l). 



Centaurus A (Cen A; NGC 5128) is the closest example of 
a giant elliptical g alaxy, located at a distance of 3.8 ±0.1 Mpc 
jHarris et al]|2010l) . It has a warped disk, likely the result of a past 
merger event, which contains an abundance of both gas and dust, 
and has a central active galactic nucleus (AGN) with radio jets 
extending out to large distances in both directions. Thus, Cen A 
is an ideal laboratory for studying variations in the ISM at high- 
resolution and to look for effects the AGN might have on the neigh- 
bouring disk. Note that while these interesting traits do make Cen A 
a somewhat peculiar elliptical galaxy, some studies show that over- 
all Cen A is still a normal elliptical on global scales, and that its 
proxi mity is what gives us insight to its peculiarity (e.g. iHarrisI 

l2oToh . 

Cen A has been observed extensively in the past over a wide 
range of wavelengths. For a comprehe nsive review o f the numerous 

S ihysic al characteristics of Cen A, see llsraej j 19981) and lMorgantil 
2O10h . Observations of H I covering a large area of the sky sur- 
round i ng Cen A (e.g van Gork om et al. 1 990; Schiminovich et al.l 
1 19941 ; iMorganti et d J 120081 ; IStruve et alj |201(I) reveal emission 
throughout the disk and in shells at large distances from the disk, 
supporting the merger scenario. The central core itself shows both 
blueshifted and r edshifted Hi absorption iMorganti et al. I l2008l ; 
IStruveetalJl2010l) , consistent with a circunmuclear disk. 

Due to its low latitude in the southern sky (13''25™27.6^ 
— 43°0r09"), only a small number of telescopes are capable of 
observing Cen A in CO, such as the Swedish-ESO Submillime- 
ter Telescope (SEST; e.g. PhiUips et al."l 987| ; lEckart et alj|l990j ; 
lOuillen et al. 1992; Rydbecketal. 1993). Detections in the J = 
1 — transition by Eckart et al.l ( Il990al) show that the emission pri- 
marily traces the optical disk, but drops off significantly beyond 
about 90 arcsec on e ither side of the c e ntre of the disk. In addition, 
lEckartetal.Ul990ah and lEckart et al. Id 1990b') report detections of 
the ^■^CO and ^^CO J — 1 — line in absorption against the nuclear 
non-thermal continuum, and suggest the same molecular clouds are 
responsible for both the emission and absorption features. 

Spitzer Space Observatory IRAC observations of Cen A 
( lOuillen et al1 l2006^ reveal a prominent parallelogram-shaped ring 
of dusty material in the central disk at higher resolution than ear- 
lier mid-infrared observations of Cen A using I SOCAM on the 
Infrared Space Observatory (ISO; iMirabel et al.|[l9 99). Observa- 
tions of Cen A at 450 and 850 p.m with SCUBA and 870 ^im 
with LABOCA show a warped 'S'-shaped structure close to the 
nucleus spanning approxi mately 5 arcmin ac r oss, similar to th e 
mid-infrared observations fceeuw etal.ll2002l ; IWeifi et al.l l2008h . 
Less-defined, fainter emission surrounds this structure and extends 
over much of the optical dust lane. To explain the 'S'-shaped disk, 
lOuillen et alj {2006) model the warped morphology of the inner 
infrared/submm disk using a 'tilted-ring' model, consi sting of a se- 
ries of concentric rings. The IRAS measurements of lEckart et al.l 
(ll990A reveal an average dust temperature of about 4 2 K, with the 
tempe rature decreasing towards the edge of the disk. iLeeuw et al.l 
( 12002 ) also find a tem perature of ~ 40 K in the disk that decreases 
outward. Furthermore, IWeifi et al.l ( I2OO8I) show the first detections 
at submm wavelengths of non-thermal emission from the radio 
lobes. Studying the thermal emission within the disk, their two- 
component dust model shows the dust has temperatures of 30 and 
20 K (warm and cold components, respectively), but no decreasing 



tre nd at larger radii. Lastly, a complementary study to this paper 
bv lAuldetai] l l2012h analyses the dust properties of the remains of 
the smaller merger component in Cen A at large radii from the cen- 
tre, u sing data from the Herschel Space Observatory dPilbratt et al.l 
[20i3). 

In spite of the number of observations of the disk of Cen A in 
both the dust and gas components, th ere are few publishe d values 
of the gas-to-dust ratio in this galaxy. IStickel et al.l j2004h report a 
ratio of ^ 300 in a region north of the disk using data taken with 
ISOPHOT on ISO, and Israel ( 1998) suggests an upper limit of 450. 
Even for other early-type galaxies, only a small number of surveys 
report gas-to-dust ratios. The SINGS survey found ratios similar 
to tha t of the Milky Way dDraine et al]|2007h , while iLeeuw et al] 
( l2008h found a range of values between ~ 230 and 400 for the in- 
ner H2 -dominated regions of seven ellipticals. In comparison th e 
typical value for the Milky Way is ~ 160 (e.g. lZubko et al.ll2004h . 
With Herschel we have the capability to study the gas-to-dust ra- 
tio for the disk of Cen A at unprecedented resolution, and we can 
investigate how this ratio varies throughout the galaxy including 
regions near the AGN, and try to understand how the AGN affects 
the structure of the ISM. Herschel also gives us increased sensitiv- 
ity compared to ground based observatories, and extends coverage 
to longer wavelengths than previous space missions, which are cru- 
cial to measuring the dust content of the ISM. 

Here we present a detailed analysis of the gas-to-dust mass 
ratio of Cen A through the central disk. We use photometric obser- 
vations wit h both the Photodetec tor Array Camera and Spectrome- 
ter (PACS; |Poglitsch et al.ll2010h and the Sp ectral and Photometric 
Imaging Receiver (SPIRE; Grif fin et al I2OI O) on Herschel, as well 
as observations with the HARP-B receiver mounted on the James 
Clerk Maxwell Telescope (JCMT). In Section|2]we present the ob- 
servations and data reduction process. In Section[3]we present our 
Herschel and JCMT photometry, our spectral energy distribution 
(SED) model and maps of the best-fitting dust temperatures, dust 
and gas masses, and the gas-to-dust ratio across the disk. In Sec- 
tion |4] we discuss the possible mechanisms that might explain the 
characteristics of our maps, and we conclude in Section|5] 



2 OBSERVATIONS 

2.1 Herschel Observations 

We have obtained PACS photometry at 70 and 160 pm (OBSIDs 
1342188855 and 1342188856), and SPIRE photometry at 250, 350 
and 500 pm (OBSID 1342188663) as part of the Very Nearby 
Galaxies Survey (VNGS; PL: Christine Wilson), which is a Her- 
schel Guaranteed Time program. All of these images cover a 37 ar- 
cmin X 37 arcmin area on the sky centred on Cen A. In Table[T]we 
summarise the basic characteristics of the original Herschel maps, 
which are shown in Figure[T] 

2.1.1 PACS Observations 

The PACS 70 and 160 /im observations were obtained simultane- 
ously in scan-map mode at the 'medium' scan speed of 20 arc- 
sec s~^. The scans were carried out in orthogonal directions to 
produce a square map with homogenous coverage, with a scan-leg 
length of 37.0 arcmin. The observations were then processed us- 
ing a pipeline that has been modified from the official one released 
with the Herschel I nteractive Processing Environment (HIPE, ver- 
sion 4.0; lot^l2010h , using calibration file set (FM,v4). Reduction 
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Table 1. A summary of the Herschel obsei'vation information as well as the pixel noise and calibration uncertainty values. The measurement uncertainty 
is pixel-dependent (see Section l2.1.3l for details) and thus values are not presented here. 



Wavelength 
(lira) 


Beam Size 
(arcsec) 


Pixel Size 
(arcsec) 


Colour Connection''' 


Pixel Noise'' 
(mJy pix^^) 


Calibration Uncertainty'^ 
(percent) 


70 


5.76 


2 


1.043 


0.009 


3 


160 


12.13 


4 


1.029 


0.060 


5 


250 


18.2 


6 


0.9836 


0.147 


7 


350 


24.5 


8 


0.9821 


0.114 


7 


500 


36.0 


12 


0.9880 


0.110 


7 



PACS values are divisive and SPIRE values are multiplicative. 

These uncertainties are for the pixels at their native pixel scale, as listed in this table. 
^ We have ignored the fact that the SPIRE calibration errors are correlated between all three bands and our total 7 percent error comprises 5 percent 
correlated error and 5 percent un-correlated eiTor 



steps include protecting the data against possible cross-talk effects 
between bolometers in the array, and deglitching only at the sec- 
ond level. This involves comparing the signal detected within all of 
the pixels in the data cube that cover the same location in the sky, 
and flagging any pixels that are outliers as glitches. In addition, 
flat-field corrections were applied to the images, and the signal was 
converted to units of Jy pix^^. We also applied correction factors 
of 1. 1 19 and 1. 174 to the 70 and 160 /im maps, respectively, to up- 
date the images with the most recent set of calibration files. Once 
these steps were done in HIPE, th e final map-ma king was done 
in SCANAMORPHO^ (version 6.0) ( iRousselllloIlh using standard 
settings and pixel sizes of 2.0 and 4.0 arcsec for the 70 and 160 /^m 
images, respectively. We applied a colour correction to each map to 
accommodate the non-monochromatic nature of the galaxy's spec- 
tral energy distribution (SED) through each of the PACS filters. To 
do so, we first measure the slope of the spectral energy distribution 
(SED) at 70 and 160 /im using our one-component blackbody fit 
(see Section lBTI for details), and applied the colour corrections for 
those slopes. Then we re-calculate the slope of the SED for the cor- 
rected fluxes and check for any change in the slope. We find slopes 
of -3 and -l-l at 70 and 160 /im, respectively, and thus we divide 
these images by correction factors of 1.043 and 1.029 (Miiller et al. 
I20T1I) . Lastly, we determine the mean sky value within numerous 
apertures overlaid onto off-target regions surrounding the galaxy at 
each wavelength, and then determine the average mean across all 
apertures. We take this mean as our background contribution and 
subtract it from each PACS image. 

The final 70 and 1 60 /im images have been convolved with 
appropriate kernels (see iBendo et alj ( 1201 2|) for details) to match 
the 36 arcsec resolution of the 500 /xm image, and then regridded 
to a 12 arcsec pixel scale, again to match the 500 /im image. 

2.1.2 SPIRE Observations 

The SPIRE observations were done in large scan-map mode with 
the nominal scan speed of 30 arcsec s^^ and the cross-scanning 
method. Data reduction was carried out in HIPE using a pipeline 
that has been altered from the standard one, in which the custom al- 
gorithm for temperature drift corrections, BRIght Galaxy ADaptive 
Element (BriGAdE), was used in lieu of the standard algorithm. 
A description of the data processing including how BriGAdE 

^ http://www2.iap.fr/users/roussel/herscheiy 



applies this correction is given in lAuld et alj Jioil) , [ Smith et al.l 
(2012) and Smith et al. (in prep). The final maps were made us- 
ing HIPE's map-making tool, and are in units Jy beam~^ with 
pixel scales of 6, 8 and 12 arcsec at 250, 350 and 500 /im, re- 
spectively. However, we converted the units from Jy beam~^ to 
Jy pix~^ for consistency with the PACS images using beam ar- 
eas for this conversion of 426, 771 an d 1626 arcsec'^ for the 250, 
350 a nd 500 /im images, respectively dSPIRE Observers' Manuall 
I2OIII) . We calibrated the flux to correct for the extended nature 
of Cen A, using multiplicative K4 correction factors of 0.9828, 
0.9834 and 0.9710 for the 250, 350 and 500 fim, respectively 
dSPIRE Observers' Manija3l2oT l|). Lastly, we applied a colour cor- 
rection to each image using the same method as for the PACS im- 
ages (see Table[T]for specific values), as well as a background sub- 
traction. The five Herschel images are presented in Figure[T]in their 
original resolutions. 

For analysis purposes, we have convolved the 250 and 350 /im 
images with a Gaussian beam to match the 36 arcsec resolution of 
the 500 /im image. The point-spread function of the SPIRE beam is 
well-represented by a Gaussian functi on down to about t hree per- 
cent of the peak in all three wavebands dGriffin et al.l2010l) , justify- 
ing our use of a Gaussian kernel. Next, these maps were regridded 
such that they all had the pixel scale of the 500 /im image, 12 arc- 
sec. 

2.1.3 Noise and Uncertainties 

There are three sources of noise and uncertainty that we take into 
account for the Herschel images. The dominant uncertainty for 
all images is the calibration uncertainty. For the PACS 70 and 
160 /im ima ges, the calibration uncertaintie s are 3 and 5 percent, 
respectively IpACS Observer's Manuall201ll) . For the 250, 350 and 
500 /i m images, the calibration uncert ainties are taken to be 7 per- 
cent ( SPIRE Observers' Manua]|201lh . There is also an underlying 
measurement uncertainty associated with the flux values in each 
pixel, and a map with these uncertainties is produced as part of 
the general reduction. Finally, there is a contribution to the error in 
each pixel from the background subtraction. To calculate this noise, 
we first find the mean value within a number of apertures over- 
laid on the convolved and regridded images, just as we did for the 
background subtraction. Then we find the average mean across all 
apertures (at each wavelength), and then determine the standard er- 
ror, given by the standard deviation of the individual aperture mean 
values divided by the square root of the number of apertures. Note 
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Figure 1. The five Herschel images at their native resolutions. The colour scale units are log(Jy arcsec" ^ ), and the beam size is in the lower left comer of each 
image. Note that the radio jets become more prominent at 350 and 500 fj,m due to synchrotron radiation. 



that we have not accounted for the confusion nois e in individual 
pixels in the SPIRE maps (see iNguven et aPbOlOh : however, we 
expect it to be trivial except perhaps in the extreme outer disk, and 
the global average of the confusion noise has been removed with 
the background subtraction. A summary of the pixel errors from 
background subtraction and calibration uncertainties is presented 
in Tabled 



2.2 JCMT Observations 

The "CO J = 3 - 2 transition at a rest frequency of 345.79 GHz 
was mapped using the HARP-B instrument on the JCMT, located 
on Mauna Kea in Hawaii. These observations were taken over the 
nights of 2010 January 24-26 as part of project M09BC05 (PI.: 
Tara Parkin), with a telescope beam size of 14.5 arcsec. We com- 
bined nine overlapping jiggle-maps (each with an integration time 
of 50 s) instead of creating a raster map due to the low eleva- 
tion of Cen A from the location of the JCMT. The maps were 



observed using beam-switching with a chop throw of 150 arcsec 
perpendicular to the major axis of Cen A. Each jiggle-map covers 
a 2 arcmin x 2 arcmin area on the sky, and our final map covers 
an area of 10 arcmin x 2 arcmin, to observe the dust lane out to 
£^25 /2, which is half the distance from the centre of the galaxy to 
where the optical magnitude falls to 25. The backend receiver, the 
Auto-Correlation Spectrometer Imaging System (ACSIS), was set 
to a bandwidth of 1 GHz with 2048 channels, giving us a resolution 
of 488 kHz or 0.43 km s'^ 

The data were red uced using the StarlinkQ software pack- 
age ICurrie et al.!'2O08'), maintained by the Joint Astronomy Cen- 
tre. We first flagged the data for any bad pixels in the raw data, and 
then built a cube (two spatial dimensions and one in frequency) 
with 7.5 arcsec pixels combining all of our observations. Next, we 
masked out detections of the ^^CO J = 3 — 2 line in our cube 

^ The Starlink package is available for download at 
|http://starhnk.jach.hawaii.edir| 
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using an automated fitting routine, then fit a third order polynomial 
to the baseline in order to remove it from the cube. We then binned 
this cube to a resolution of 20 km s~^ to measure the average rms 
in the line free regions, which we determined to be = 16 mK. 

The first three moment maps (integrated intensity, integrated 
velocity and velocity dispersion) were created using the routine 
'findclump' (part of Starlink), which looks for detections of 
greater than 2a. We first convolve the unbinned cube with a Gaus- 
sian beam to match the 36 arcsec resolution of the Herschel SPIRE 
500 fim map. Next, we determine the average rms within the line 
free regions of the spectrum and expand the two-dimensional noise 
map into a cube. Then we create a signal-to-noise cube by divid- 
ing the convolved cube by the noise cube and apply 'findclump' 
to it. This method allows us to make detections of low-level emis- 
sion better than just applying 'findclump' to our image cube. The 
routine creates a mask of the detections in a smoothed version of 
the cube and then multiplies this mask by the signal-to-noise cube 
to obtain a cube containing only detections above several differ- 
ent a cutoffs. Then this cube is collapsed along the velocity axis 
to create the three moment maps in terms of the signal-to-noise. 
Finally, the integrated intensity-to-noise map is multiplied by the 
two-dimensional noise map to obtain the true integrated intensity 
map. For a more de tailed description of our reduction method see 
IWarrenetal.l ( l2O10l) . 

To evaluate the uncertainty in the integrated intensity map, we 
make use of the following equation: 



^I=[^voVN^Jl + P^, (1) 

V "base 

where Aw is the width of each channel in km s^^, a is the noise 
in units of Kelvin, A'line is the total number of channels within 
the spectral line itself, and A'^basc is the number of channels used 
for fitting the baseline. Several calibration sources measured during 
these observations were compared to standard spectra, and from 
these we estimate the calibration uncertainty to be 10 percent. The 
integrated intensity map is presented as contours in Figure|2] 

2.3 Ancillary Data 

The H I map of Cen A was kindly provided to us by Tom Oosterloo 
and previously published in Struve et al. (2010). This map has a 
resolution of 19 arcsec and units of 10^^ cm~^, with a pixel scale of 
4 arcsec. We smoothed the image to a resolution of 36 arcsec using 
a Gaussian beam to match the resolution of the SPIRE 500 
image, and regridded the image to match the 12 arcsec pixel size. 

We also have a map of the radio continuum at 1 .425 GHz, 
originally published as part of the IRAS Bright Galaxy Atlas 
dCondon et alj|l99^ . and retrieved from the NASA Extragalactic 
Database (NED). This map allows us to trace the jets extending 
out from the central AGN. We convolved this radio image with a 
Gaussian beam and then regridded it to match the SPIRE 500 
image. 



3 RESULTS 
3.1 SED Fitting 

To determine the dust temperatures within the disk of Cen A, we fit 
the far-infrared and submm part of the spectral energy distribution 
(SED; 70 - 500 /im) with a simple modified blackbody with P — 2, 

I{u,T) = C xu^ X B{u,T), (2) 



where C is a scaling constant to match the model to our observed 
fluxes. The parameter /? originates from the dust emissivity func- 
tion, Ki, — Ko[v / Uq)^ . 

To constrain our modified blackbody fits, we created uncer- 
tainty maps for each wavelength that are a combination of the mea- 
surement uncertainty in the flux value of each pixel, the background 
noise in each pixel, and the calibration uncertainties for both PACS 
and SPIRE (see Section |2. 1.3| for details), added in quadrature. The 
value in each pixel of the resulting map is taken to be the total un- 
certainty for the equivalent pixel at a given wavelength. 

For each pixel in the images that has a S/N 10 at all five 
Herschel wavelengths, we fit Equation[2]to our flux measurements 
to determine the best-fitting temperature T and constant C for that 
pixel. An example of a SED fit for a typical pixel is shown in Fig- 
ure|3] and the pixel location is shown with a cross in Figure|4] We 
also show the global SED that was obtained by summing the flux at 
each wavelength within all of the pixels with good fits in Figure[3] 
and fitting the resulting totals. Note that we tested the effect of re- 
moving the 70 /im data point from our SED fit to ensure it was not 
forcing the curve to peak at shorter wavelengths, and that it was not 
from a separate thermal c omponent of the dust ISmith et al.ll20ld: 
iBendo et al.|[2010all2012h . We found that the 70 ^m flux falls on 
the best-fit curves within uncertainties for these tests and thus we 
use all five wavelengths in our fitting routine. We then create a map 
of the dust temperature where pixels with poor chi-squared fits at 
the 95 percent confidence level have been masked out. In general 
the pixels with poor fits coincide with those where non-thermal ra- 
diation is making a significant contribution to the flux at 500 /im. 
The temperature map is presented in Figure |4] with 500 /^m con- 
tours overlaid to point out the location of the non-thermal emis- 
sion. The temperature is about 30 K near the central AGN region, 
and then smoothly falls off to about 20 K in the outer disk, and the 
average temperature throughout the entire region is approximately 
(21.8 ±0.3) K. 

We have also tried allowing both /3 and temperature to vary 
in Equation |2] to investigate if the model fits to our data improve. 
While there is some small variation in jS throughout the disk, these 
deviations are small and the average value of /3 is 2.07 ± 0.07. 
Furthermore, there are no obvious trends in the variation of /3 with 
radius. Thus, the results of our test are consistent with our one pa- 
rameter fit, giving us confidence in our fixed /? fits. 

We evaluated the uncertainty on the dust temperature by im- 
plementing a monte carlo algorithm on our modified blackbody 
fitting routine. For each wavelength, an artificial flux value lying 
within the uncertainty bounds of our observed measured flux value 
is generated randomly using a Gaussian number generator. Then we 
run the modified blackbody fitting routine on the simulated dataset 
and extract the temperature and constant for the best fitting function 
corresponding to that dataset. This process of creating artificial data 
and then fitting a function to them to determine the temperature is 
repeated 800 times for each pixel. We take the standard deviation of 
the temperature distribution for each pixel to be the uncertainty in 
the temperature. This uncertainty varies throughout the map, with 
an average of about I percent per pixel. The highest uncertainties 
are 5 percent around the very edge of our map, where the signal to 
noise is lower for the flux maps, and we note that these uncertain- 
ties exclude the uncertainties arising from our assumed value for /3 
and thus emissivity. 
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Figure 2. Hi map of the disk region of Cen A in units of 10^^ cm~^ from lStruve et alj i201(]|l (colourscale), with a beam of 19 arcsec. JCMT CO J = 3 — 2 
contours with a resolution of 14.5 arcsec are overlaid in black. CO contour levels are 0.4, 0.8, 1.6, 3.2 K km s^^. . . 25.6 and 38.4 K km s^^, where the 
temperature units are TpJ*. 
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Figure 3. left: The SED for a typical pixel (pixel location shown as a cross in Figure|4). The best-fitting modified blackbody is shown as a solid line, while our 
measured flux values are represented by the data points, right: A global SED for the entire disk, consisting of all pixels contributing to the total dust mass (see 
Sectionim. 



3.2 Dust Mass 

The equation we use to evaluate the dust mass at a specific fre- 
quency is 



Afdust 



(3) 



when the monochromatic flux is optically thin. In this expression, 
Sv is the flux from the source, D is the distance to the source, 
Biy, T) is the equation for a blackbody and is the dust emissiv- 
ity. Here we adopt a value for the dust e missivity at 25 pm, K250 , 
of 3.98 cm^ g"^ from the dust model of lDraind j2003lfl . Substitut- 
ing ^250 into Equation[3l the dust mass becomes 



q / 57.58 K \ 

5.258 X 10-^ 5*250 (e^^ - l) Mq. 



(4) 



We checked that the SED of Cen A is optically thin by assuming 
that the pixel with the highest flux at 70 pim is a worst-case scenario 



^ These data are also available at 

|http://www.astro.princeton.edu/~draine/dust/dustmix.html| 



in terms of optical depth r, and measured the flux in the same pixel 
at 250 /im. Using the results from our SED fitting we calculate 
the dust mass within that pixel using Equation [4] and convert the 
dust mass to a dust mass surface density E using the area of one 
12 arcsec pixel. The optical depth at 250 /im can be calculated as 
r25o = E«;25o. Using this method we find that r25o ^ 1, thus we 
can assume the galaxy is optically thin at 250 /im. This method can 
be repeated for the other wavelengths as well. 



Substituting our temperature and the SPIRE 250 ^va fluxes 
returned by our model SED fits into Equation|4l we obtain a pixel- 
by-pixel map of the dust mass, which is shown in Figure |5] The 
dust mass distribution peaks toward the centre, and falls off in the 
outer regions. Uncertainties range from about 3 to 30 percent with 
the lower uncertainties corresponding to pixels closer to the centre 
of the disk. The typical uncertainty in each pixel is ~ 5 percent. 
Summing up the individual dust masses within each pixel leads to 
a total dust mass in the disk of (1.59 ± 0.05) x lO'^ Mq. As a 
check, we also evaluated the total dust mass using the global SED 
fit results, and obtain a total of 1.47 x 10'' Mq. 
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Figure 4. Temperature map of the dust disk of Cen A (colour scale), in units of degrees Kelvin. Contours of the 500 /im SPIRE map are shown in grey for 
reference. The contours for the 500 ^m image are 0.0139, 0.0347, 0.0694, 0.347, 0.694, 1.39, 2.78, 4.17 and 4.86 mJy arcsec"^ -phe black cross is the location 
of the pixel for which the best-fitting SED is shown in Figure|3] Note that the central region does not have any reasonable fits due to the non-thermal emission 
present in the 350 and 500 ^m maps. 
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Figure 5. Dust mass map of Centaurus A. Units of the colour scale are M© pc . Contours of the 1 .425 GHz radio continuum are overlaid in black to show 
the regions where synchrotron radiation is present. 



3.3 Gas Mass 



XcO^CO{3-2) / -^CO{3-2) 



77mb 



'CO{1-0) 



(5) 



We convert our CO J = 3 — 2 integrated intensity map to a column 
density of molecular hydrogen using the following equation: 



where Xco is the C0( J = 1 — 0)-H2 conversion factor, /co(3-2) 
is the integrated intensity of the CO J = 3— 2 emission, J7mb = 0.6 
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is the conversion factor from an antenna temperature Ta* to a 
main beam temperature Tmb for the JCMT, and the factor within 
the brackets is the ratio of the CO J = 3 — 2 line intensity 
to the CO J — 1 — line intensity. Despite the additional un- 
certainty introduced by this line intensity ratio, we choose to use 
our CO J = 3 — 2 map in this analysis because of its better 
intrinsic resolutio n (14.5 arcsec) and s ensitivity compared to the 
CO J = 1 - of lEckart et ai] h990ah . We assume a Xco factor 
of (2± 1) X 10^° cm~^ (K km s'^y'^, typical for the Milky Way 
( IStrongetalJl988h . and a CO J = 3-2/CO J = 1 - ratio of 0.3, 
which is a suitable ratio for the diffuse IS M as found by the J CMT 
Nearby Gala xies Legacy Survey (NGLS dWilson et alj|2009h . We 
note here that i Wilson etal] ( |2009|) and other groups from the NGLS 
have seen ratio variations within other galaxies by up to a factor of 
2-3, which may increase our molecular gas mass uncertainties (see 
Section|4]for further discussion). 

Next, we convert our column density to a molecular gas mass 
map via 

= ^pixA^Ha^TlHa (6) 

with ^pix the area of one pixel (12 arcsec x 12 arcsec) at the dis- 
tance of Cen A, and the mass of one molecular hydrogen 
atom. The total mass of molecular hydrogen, over the entire cover- 
age of our map, is (1.42±0.15) x 10^ Mq, where we have excluded 
uncertainty contributions from the assumed CO J — 3 — 2/CO J = 
1-0 ratio. 

Following the same method we also convert the units of the Hi 
map (Figure O to units of mass. Then the molecular gas mass and 
atomic gas mass maps are combined in such a way that in regions 
where we have coverage for both Hi and H2, the total gas mass is a 
sum of the two, and for the remaining regions the mass is from Hi 
or H2 alone. We then multiply the resulting map by 1.36 to obtain 
the total gas mass, including helium. Summing over all pixels, we 
find a total gas mass of (2. 7± 0.2) X 10^ Mq. Our map of the total 
gas is presented in Figure (6] and covers about half of the area for 
which we have dust mass measurements. 

3.4 The Gas-to-Dust Mass Ratio 

We present the gas-to-dust mass ratio in Figure|7] which covers all 
pixels for which we have both a gas mass and a dust mass. The box 
outlined in black shows the region that has coverage both in Hi and 
H2, which is 10 arcmin x 2 arcmin in area. As is seen in this figure, 
the majority of the disk has a gas-to-dust ratio similar to that of the 
Galaxy with an average value of 103 ± 8; the Galactic value varies 
from a bout 120 (,Li & Dra ine 2001) to ~ 160 (Zubko et al. 2004) to 
~180 jPraine et al J 20071 and references therein). Interestingly, the 
gas-to-dust ratio increases to a significantly higher value of almost 
300 in regions closest to the AGN. Typical uncertainties in this ratio 
are about 10 percent. We also present the average gas-to-dust ratio 
as a function of distance from the centre of Cen A along a position 
angle of 30 deg north of west in Figure[8] This plot emphasizes the 
radial trend of the gas-to-dust ratio along the disk. In other galaxies 
such as those from the SINGS survey, th e gas-to-dust ratio varies 
from ~ 136 to ~ 453 jPraine et alj2007h . We discuss the possible 
origin of this high central gas-to-dust mass ratio below. 



4 DISCUSSION 

Our total dust mass of (1.59 ± 0.05) X 10^ Mq is comparable 
to other published numbers for Cen A in the literature; however. 



i.; ' 'iiii 
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Figure 8. A radial plot of the average gas-to-dust mass ratio along the tnajor 
axis of Cen A with a position angle of 30 deg north of west. We measured 
the average gas-to-dust ratio value within 20 bins along this position-angle, 
each measuring 2 arcmin X 0.5 arcmin in size. The red data points represent 
the bins in which we had at most 10 good pixels out of a possible 25, and 
therefore are not as well sampled as those shown in black. 



we note that this value is likely a lower limit as regions where 
non-thermal emission contaminates the dust continuum had poor 
model fits. Several publications have used two-component modi- 
fied blackbody modelling combined with far-infrared and submm 
observations to derive a dust mass for Cen A. Leeuw et al. ( 200^ 
used a wavelength range of 60-850 fim for their SED fit and found 
a total mass of 2.6 X lO'* M0 (adjusting for our distance) within 
an elliptical region 120 arcsec x 450 arcsec in size, but excluding 
the unresolved core at the centre. Since this ellipse covers a similar 
area to our observations , the discrepancy is likely due to an inade- 
quate chop throw in the lLeeuw et al.l j2002h 850 dataset, lead- 
ing to difficulties quantifiying the background. Our observations 
show faint emission at 500 /im extending beyond the 120 arcsec 
chop throw the authors used, and oversubtracting the background 
will lead to a lower dust mass estimate. WeiB et al. (2008) studied 
the submm continuum using LABOCA 870 fim data and combin- 
ing these data with fluxes spanning 25-870 fim they found a higher 
mass of 2.2 x 10^ Mq. They subtract the flux contribution of the 
unresolved core from the total flux prior to estimating the dust, but 
do caution that they may not have removed all of the non-thermal 
emission arising from the radio jet. 

The gas mass we have obtained from our H2 map alone is 
(1.42 ± 0.15) X 10^ M0, incl uding the central reg ion which we 
masked out for our dust fitting. Eckart et alj ( Il990al) determined a 
molecular gas mass of 1.9 x 10* M0, after we adjust t heir result 
for ou r assumed Xco factor and distance. In addition, iMorgantil 
( l2010h suggests th e total molecular gas mass is 4 x 10* Mq base d 
on observations of lEckart etaEI jl990ah and lOuillen et alj ( Il992h . 
Our molecular gas mass is thus higher than those values previously 
published; this is likely due to the improved sensitivity of our ob- 
servations. The m olecular gas mass i s about twice the atomic mass 
of 4.9 X 10* Mq ( IStruveetal.ll20 1Q^. 

The high gas-to-dust mass ratio towards the central AGN is 
a most interesting result. Several studies of spiral galaxies have 
found that the gas-to-dust ratio is not constant throughout the disk, 
but rather decre ases with decreasing distance from the centre of 
the galaxy (e.g. Muiioz-Mateos et alj |2009| ; iBendo et alj l2010bl; 
Magrini et al...201 1.) . in contrast to what we have found here. We 
have made a rough estimate of the CO J = 3 — 2/CO J = 1 — 
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Figure 6. The total gas mass distribution of the disk of Cen A. This map includes both the Hi mass distribution, as well as the H2 mass, where we have 
coverage for it, and has been multiplied by 1.36 to account for helium as well. Units for the colour scale are Mq pc~^. 
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Figure 7. Gas-to-dust mass ratio of Centaurus A (colourscale), with PACS 70 /xm contours (at original resolution) overlaid in black. Contour levels are 0.75, 
1.5, 2.5, 7.5, 20, and 37.5 mJy arcsec"^. The box shown in black highlights the region in the disk for which the gas contribution is from Hi and H2, and is 
10 arcmin X 2 arcmin in size. 



ratio throughout the central disk using the CO J = 1 — map of 
lEckartetal.l ( ll990iJ) to check if our assumption of 0.3 for this ra- 
tio is valid in the centremost regions of the disk. Our calculations 
show an average value of 0.35, and more importantly, they show no 
significant variation with radius. Furthermore, looking at Figure [7] 
we see that there is a smooth transition of the gas-to-dust mass ra- 
tio between the H2-dominated region to the Hl-dominated region in 
the outer disk, which also suggests CO excitation is not the primary 
cause of the high gas-to-dust ratio. 

The Xco factor is one variable that might affect the gas-to- 
dust ratio if it is not constant throughout the disk. We have as- 
sumed the Galactic value for our calculations; however, studies 
have shown th at this conversion fa ctor can vary w ith metallicity 
(IWilsorJll99l: Jsrael 1997, 2000: Barone et aLlfioO Q : Israel et al. 
I2OO3I: IStrong et alj |2004 llsraell l2005l: iLerov et al.l I2OI lb If the 



metallicity decreases as a function of radius within Cen A, we 
would require a larger conversion factor at larger radii, which in 
turn would increase the gas-to-dust ratio farther out in the disk, 
potentially bringing it up to a similar value as that in the cen- 
tre. Alternatively, the Xco factor ca n be smaller by up to a fac - 
tor of ~ 4 for starbursting galaxies jPownes & Solomorj[l998h . 
This effect could reduce the gas-to-dust ratio in regions of strong 
star forming activity; however, the total far-infrared luminosity 
for Cen A is ^ 10^° L©, while the Xco conversion factor is 
typically affected at luminosities greater than ^ 10^^ Lq (e.g. 
iDownes. Solomon & Radfordlll993l ; rSolomon et aljl997h . 

One possibility is that dust is being destroyed by the jets ex- 
tending out in either direction from the AGN. It is plausible that 
some of the dust grains are being destroyed by the AGN either 
through dust sputtering or through shocks. The total X-ray lumi- 
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nosity of the AGN of Cen A is 8 x lO"^ erg s"^ jRothschild et all 
l201lh , w hich is comparable in X-ray luminosity to the AGN found 
in M87 iWemer et al]|2006h . It might be possible for these pho- 
tons to d eplete the dust in the surrounding ISM. In support of this 
scenario. iBaes et J j2010 ) did not detect thermal dust emission in 
M87, which also contains a jet. 

It is also possib le that the dust is entrained in the jets of Cen A. 
iRoussel et al] bOld) found a dusty halo surrounding the starburst- 
ing galaxy M82 using Herschel observations. They considered en- 
trainment by a galactic wind as a possible explanation even though 
they finally concluded that the origin was likely due to tidal in- 
teractions. We cannot rule out entrainment by galactic winds as a 
possibility for Cen A. Unfortunately, we cannot probe close enough 
to the jets to further investigate this assumption due to our limited 
resolution and the non-thermal emission component. 

Alternatively, a larger gas reservoir in the centre could also 
increase the gas-to-dust ratio. This would require dust-poor gas 
falling onto the disk and being funneled toward the centre. How- 
ever, if this was the case, we would likely have found a high gas-to- 
dust throughout the disk and not concentrated toward the innermost 
part of the disk. It also does not seem plausible for the gas to be mi- 
grating through the disk faster than the dust. 

Interestingly, there is a large ring of material that is seen in the 
mid- and far-infrared close to the regions in which we see a high 
gas-to-dust ratio. MirabeLelLal- 1 1999) first noticed the 'S' shaped 
structure with ISOCAM observations at 7 and 15 /im, and postu- 
lated that this structure was a bar at the centre of the galaxy. Later, 
observations jLeeuw et aLll2002l : lOuillen et al.ll200d) supported a 
tilted ring scenario to explain the structure observed. Our Herschel 
observations at 70 ^m also show emission from this ring (Figure[Tll, 
and in fact, the high gas-to-dust ratio appears to correspond with the 
70 ^,m contou rs of the ring, as shown in Figure [7] The models by 
lOuillen et al] (^006) suggest that there is a deficit in the dust dis- 
tribution at radii less than about 50 arcsec (~ 920 pc). This is at a 
similar radius to our observed high gas-to-dust ratio. A dust deficit 
throughout r ^ 50 arcsec might seem less likely to be related to 
the AGN and is perhaps better related to star formation activity. 
Observations probing closer to the nucleus are necessary to deter- 
mine if the gas-to-dust ratio remains high or increases further as we 
approach the AGN. 



5 CONCLUSIONS 

We have presented new observations of the dust continuum from 
the PACS and SPIRE instruments on board the Herchel Space Ob- 
servatory at 70, 160, 250, 350 and 500 ^m. In addition, we present 
new data tracing the CO J = 3 — 2 transition taken with the HARP- 
B instrument mounted on the JCMT. We have used these data to 
probe the interstellar medium within the disk of Centaurus A on 
a pixel-by-pixel basis. We observe the ring of emission at 70 /xm 
previously reported at infrared and submm wavelengths, while our 
500 pm (and to some extent the 350 ^m) images show detections of 
th e non-thermal cont inuum emission previously reported at 870 
bv lWeiB et al] ( |2008|) . 

We model the dust spectral energy distribution using a single 
component modified blackbody, and find temperatures of ^ 30 K 
in the centre that decrease with radius. Using the temperature map, 
we find the total dust mass to be (1.59±0.05) x 10^ M©. The total 
gas mass is (2.7 ± 0.2) x 10^ M©, and combining the dust and gas 
masses we have produced a gas-to-dust ratio map. The average gas- 
to-dust ratio is approximately 103 ± 8, similar to that of the Milky 



Way, with an interesting peak to about 275 towards the centre of 
the galaxy. After exploring several possible scenarios to explain 
the high gas-to-dust mass ratio, the most appealing one is that of 
a correspondance between the high gas-to-dust mass ratio and the 
ring about 1 kpc in size , whic h is best fit by a warped tilted ring 
model by Ouil len et al. I ( l2006h that consists of a deficit of dusty 
material inwards of the ring. The deficit of dust may be due to X- 
ray emission from the central AGN destroying the dust grains or 
the jets are removing dust from the centre of the galaxy. 
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